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ABSTRACT 

LP emserCUu- memati more arises from a Poisson distribution 
with parameter A, when the distribution of A itself is of 
the form BN wea) ee Acel[0,b], where c is a scaling factor 
Miceomen are Surieuly POSitive parameters. However, the 
mirertonal form ef Mee resulting unconditional distribution 
memos Darvtviculartyetractable, hence the study of the sta- 
meovrcal properties» ol the unconditional distribution is 
imemaced in the study of its mean and its variance. 

In regards to the modelling of a real situation, an 
estimation procedure of the parameters involved in 
Pe (p-1) 7 is discussed and a closed form of the proba- 
Pality GiIsterioum@men is derived. In addition, when accuracy 
is desired a numerical analysis of the probability distribu- 
tion is also presented. The development of the results is 
continued in Appendix A, as a preparation in computerizing 
MMemcaloulat i OM. 

Prieeeeeran application to real data is discussed for 


bac Purpese Of illustrating the model. 
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eee SRODUCTION 


Let X be a Poisson variate with the probability distri- 
Pucpston depending on a real parameter A, i.e. 
ae 


(aay) Pr{X=n|A=)} = rt = Onlin? Sia. 
Tae 





where A itself is a positive random variable assumed to have 
Peeererer Gistribution specified by the probability density 


Gemover a range of definition R. 


21s) Ses 
g(A) 


6) otherwise. 


f e@(A) dA = 1. 
R 


Depending on the functional form of g(A), the uncondi- 
ier Oria I probability Gast riouvion Of XxX may result in a wide 
variety of forms. 

With regards to the application the use of (1.1) in the 
field of accident statistics, and acceptance sampling plans 
for manufactured articles based on the counting of defects, 
Sverre cewell Knowi. The variations in the accident 
liability from individual to individual in the former case; 
and in the latter an inevitable and continuous change in 
Monier condst10ns, lead to fluctuations in the Poisson 


parameter involved. 





However, the true behavior of this random variable is 
Seecnown. Thvs, che distribution of A used so far are 
meghane more CtChamweapproximations. 

Piece cemencnn > 15 PeasOnable Co estimate the 
mersc> theec Or four moments of the distribution to get some 
general idea. 

Referring to accident statistics, this method was ex- 
plained and employed by Newbold in 1927 [1]. The results 


byemer data stated that the distribution of A is 


(k/m)* k-1_.=(k/m)A 
le ) 


G2) g(r) = he [05@) 
where k, m > Q. 

That is a Pearson type III distribution as originally 
suggested by Greenwood in 1919 [1], which is also known as 
the Gamma distribution with parameters k and k/m. 


igeuomlows that the unconditional probability of X equal 


eon is: 
Pr{X=n} = f Pr{X=n|A=A}g(d) da 
0 
k co 
Gp) _ (k/my*  Y ntk-1,-(k/m 41) gy 
ae kel 20 
Tin (ey () 
mirck mtk mtk 





The random variable X is then Negative binomial distri- 
buted with parameters k and m/(m + k). 

Using the method of moments, to estimate these parameters, 
and let k and a muomrne oc tamates Of k and m: x is the esti- 
mate of the mean, 46 is the estimate of the variance of X, 


we have 


> 


i, o se 


Bele = x) 


nN > 
iH 


In regards of modelling a true situation, the above 
fie@etewith A given A Poisson distributed and prior distribu- 
tion Gamma (k,k/m) is known as the Negative Binomial Model 
imacciGent Statistics. 

hipwovwine tmis model, many workers in the accident sta- 
tistics field achieve a good fit to real data. As an example, 
iepevecmeoncuder the following data, referred to accidents 


met at work by 122 shunters during a period of 6 years [1]. 


xX Observed Frequency 
0 40 
il 39 
26 
3 8 
4 6 
5 2 
6 i 





We now assume X is Poisson distributed with parameter A 
and Ais distributed Gamma (k,k/m), then apply the Negative 
Pamemial model touche data. 

Using the method of moments, with x = 1.27 and 5° = L642 


to estimate k and m, we have 


> 


meo= Loe 7 


m™ 2 
i 


4. 36 


hence k/m —aTecpememceonen Gamma (41.360, 3:43) distributed. 


Computing the Pr {X=n} by (1.3), and the expected fre- 


m@iencaes for n = 0,1....,6 results in 

X Pr{X=n} Expected Frequency 
0 58 2 hs SBS 
1 eZ Ti BS) 0 5H 
eS De Sys 230 1 
S: -0933 i she 
4 SOS 4.72 
5 0146 ec 

> 6 0077 93 


Performing a chi-squares goodness of fit test, results in 
the value of xo as 1.589 with 4 SeROOE of freedom, hence 
p = 80%. 

However, the Negative Binomial model does not always 


give satisfactory results as in the above example. Since 





the density function specified by (1.2) is simply an approxi- 
oeret@mecmerNenoUL UNKNOWN Situation. in most of the cases, 
the approximations are based on the facts that the random 
Variable is positive [2], and that the exponential term is 
included, a simplification is expected in subsequent compu- 
tations. However, the latter factor also makes the assuming 
iiesemiouciOon rigid in the shape. Having two parameters, 
but the shape of the Gamma distribution is mostly dependent 
on one and can only graph three typical curves as illustrated 
moan le) |. 

Thus, (1.2) may not reflect well the true behavior of 
the case where the random variable is distributed in a 
Gitierent manner, such as uniformly. 

Mockinewemecullis assumptvion, the purpose of this thesis 
femomraccvemptemvo look for a more flexible prior distribution 
of the parameter A, then study in detail the resulting 


mmeenmeativional distributiongofl the variate xX. 





TABLE I 


Typical GAMMA curves 
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Wier norOokD PRIOR DISTRIBUTION OF ) 


Define a density function g(d) of the random variable A 


as follows: 


2.1) 


a eo me ee LO. Os 
g(a) = 


O 3; otherwise 


miememo,s and b are Strictly positive, and c is a scaling 
f/rever. (i) defined as in (2.1), is assumed to have the 
fom OL the Bela distribution, but the domain of definition 
momO,b, instead of [0,1]}. For this reason, from now on 


A is said to be Beta distributed. 


eens, ona PE OF THE DISTREBUTION 
Let label (B) be the curve representative of g(iA), and 


consider the derivative of g(A) with respect to i. 





(Ce) = PeMEeee GEc(uey macs = (ate—2) 1}. 
d _ s a-l 
(25) Ap wee Cee = D5 


Depending on a and 8, Nes eammveuwi ulm Or CuUbsidae of [0,b i]: 
If LEO sb J then g(\) has a maximum (global) or minimum 
Comopoieimem [Osbls. If i [0,b] then g(A) is monotonically 


increasing or decreasing in [0,b]. 


ial 





In general, the shape of (B) near the origin is determined 
by the factor een 

Take the limit of (2.1) and (2.2), when A approaches 
zero from the right. It can be seen that (B) will be tangent 
mome@nemver pical axis for a < ior anvereeny, bLnas axis Lor 
oes, symmetry, the factor Coo = determines the shape 
eiecB)) wear b. 

Simdlarly, when A approaches b, (B) will be tangent for 
Meeeeeer anverceplL tor 6B > 1 the vertical having the equation 
X = b. 

Thus, for different combinations of ao and g, (B) may 


graph in a wide variety of curves. Some of the typical ones 


are shown in Table II, Table III, and Table IV. 


Pee olATISOTICAL PROPERTIES OF THE DISTRIBUTION 


By definition of a density function, g(A) must satisfy 


b 
jermenen). AGhu= 1. 
0 


or in this case 


ie 
c fs x°-t (p-a) FI Slie= 
0 
Thus, 
ae if oiels) 
ptte-t aca) r(g) 
and 


12 





PAB fT 


iyescal SETA curves for a < 1 





1.3 





TAD Seal 


tyomeemmaslA curves for a = 1 





14 





TABLE IV 


iyobtedimoelnaecurves for a > 1 





ine 





T(atBp) Cia Bel 
e(A) = + 1 (pe) 
Doe ra) T (a) 


Let EC A‘} be the kth moment about the origin of A, 


b 
epa*y e TCatsB) f ,otk-1 


B-1 
= (b-A) ov,\ 
pb” © “Tajr(s) 0 


: pk TC atB) IP Catk) 
>T(a) I CatBtk) 


mr perticular, 


a TCatB)T Cat) 
(2.4) ELA] = b FCayFCatBF1) 


and 


eee. 2 SCOR) SG. ) 

aoe ELA] = Bo Pay rCatete 
| _ be oy Uemeilyy, 

(ate ) atBtl) 


ipeacimbewo mye 2s) sand by (2.5) that 


aCatl) 2 OL 


2 
(2.6) ae cmomicrae ly) | CO 
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On the other hand, for all te(-ty sto) the moment genera- 


ting function of the random variable A is 


T (at® ) 


b 
pe i ae Ge 
Me Ua) BCR) 0 


A 
(267) M,(%) = Ele**] = 
mien 1S known to exist but not particularly tractable [3]. 
ieee Uncondici1onel Distribution of X 
Multiplying Pr{X=n|A=A} by g(A) and integrating over 


me}, the unconditional probability of X equal ton is 


given by 
b ,atn-1l B-1 
Pr{X=n} = rere f Eerie = Ee. dr 

Db hea Te) 10 n! 

= cre) 

7 a+B—-1 EMC, 
n! b nea Ce) 

Pale tcineg tv = "=|. ana by substituting a by (atn) in 


fe mand mocing that the functional forms of Pr{xX=n! and 
(2.7) are now similar. It is then immediate from (2.7) and 
the argument following that Pr{X=n} is not tractable. 

However, the existence of Pr{X=n} can be alternatively 
proved as follows. 


Take the sum of Pr{X=n} for n = 0,1,..., which results in 


4 Pr{X=n } = q ae ae Ps ease oe ae aan GA 
n=0 n=0 b°*B-lrrayreg) oo as 
= Tats) ya-dey yy B-1g-Ay A gy 
~ patp-l - ea 
b Heel te) 10 n=0 


ey 





By emOoting that 


00 n 
y aes = er 
ete av 
so that 
a T(atB) : a= 1 B-1 
XY Pr{X=n} = TEs a: (b-A) dA 
ae. b Tr(a )TCB) oO 
ce. ) = Pr{id<b} = 1 


ia addition, 


Iv 


0 Ae LO,b] 


eiobying that 1(a,8,n) is greater or equal to zero. Then it 
moelows DY (2.0) that Pr{xX=n} exists as a finite positive 
GuameLty . 
A closed formula of PriX=n} will be derived later in 
mec v10n LIL. 
2. Estimation Of The Parameters 
Suppose observation is made on a heterogeneous popu- 
lation (S) a mixture of t homogeneous sub-populations during 
a certain period. The observed values of X are recorded, 
and denote Xv the maximum of these values. 
Let ds be the value of i of the ith sub-population, 


and let 


18 





An = max tr, } oa ae Me ae aaa ee 
al 
and (M) be the corresponding sub-population. 
Clearly, by the expression of g(A), b is identically 
Pomou now ce unat sim epractice 


egual to pd However, A 


M° M 
it is estimated as follows. 

Assume for a moment that the members of (M) can be 
identified from the heterogeneous population (S). Let m 


be the size of (M) and Y be the Poisson variate X, particu-~ 


larly referring to (M), it is immediate that 


ee =ic)) ELY] = Avy 


Now let 5 for i = 1,2,...,m be the observed values 


Seer corresponding to the ith individual of (M). Noting 


that the individuals in (S), with the observed value X,, very 


M 
likely belong to (M) then we can assume 


(2.9) ye = Kk , for some i; 


Using the method of moment to estimate A results in 


M? 


Further, based on the assumption of continuity made 


on the random variable A, we may repartition (S) into ttl 


Wee 





homogeneous sub-populations by splitting (M) into two sub- 


Mepulations (M1) and (M2) in the manner so that (M2) contains 


a unique gun member of (M) having Y, = Xv Estimating A 


in (M1) and (M2), we have 


ys See 
lq = =p 8 


izj + 


nw 


maus, A and more general 


m2 7 Amy 


Mine 22 De Orn ees: Da ge M—I MEL et. 


“A 


Hence, Avo can be used as an estimate of i 


the estimate of b is 


Me? hee tollows that 


0) b= X 


Now let E[X] be the mean, and V[X] be the variance 


of X. By elementary reasoning often used by Bayesians; 


ECX] = ECE(X[A)] 
= E[A] 
V[X] = ELV(X]A)] + V[E(X|A)] 


= ECLA] + VCA] 


20 





Mieus, by (2.4) and (2.5) 


Hy 
ox 


ELX] 


and 


2 


aCatl ) 2 OL tp 


Le 
VEX] = b Taw yatB+1) ~ ee atp 


Denote x, s°, a and 8 as the estimates of E[X], 


V[LX], a and g respectively, then 


ay 


wN 
<= pe | and 








earls 
AN 2 ‘ AN oe. ray 
e251) 3° = Roweatatl) - oS ee fp 
CatB) Catgt+l1) (atB) atBp 
2 


Supearename x from (2.11) and dividine by x 


“ww ay mn 2 = 
@ 12) (atl) Cats) -~l = S — A 
aCatgt+l1) a 
However, 
(2-13) ars = pk 
el x 


Pcie oy esmesetctubing (2.13) into (2.18) and adding 1, 


2. 2 a becenes 


el 





(2.14) bCatl1) 2 ae -~ x + x° 
x(atBt+1) xc 


Now, from (2.13) 


OG tobe oO A 
xX 
or 
(2.15) eee =e OS 
xX 


Mae by substituting (2.15) into (2.14) to get 


Ei6) peo <A 
lOfel ae oe x 
or 
Cee) x° = (bo Fae (some x + x°) 
Selvineg fo0r @ in the above equation, it comes out 
em 5 
poe Cnet X= eee Ss) 
rAd) 0 pe rN E as 
b(s”- = X) 
SMbewituvane the value of a into (2.11), to have 
(2.18) g = a(—- 1) 


IS 


ae 





Now observe that estimating by by X implies 


M? 


(2.19) baa or 


> | [O° > 
V 
+ 


Furthermore, considering (2.9), this equation is symbolically 


written as 


eee VEN) + Bf x | 


Noting that the variance of a random variable is always 


positive, then 


VLX] - ELX] = VLA] > 0, or equivalently to 


aa x Le 
2. 20) a6 = xX = bo IC: ae - b° as a 0 
(atB) (a+B+1) (at+B ) 


Micomonsernve that by (2.17), (2210) and (2.19), if 
Geis iepative (positive) it follows that gp is negative 
(positive) and vice-versa. Suppose now that a and g both 


be negative. This implies that 


a +l o} 
°C  —_—_—_s. 
ao aop tok a + B 
since 
a ae oe 
cour 6 


23 





and leads to 


aA wn oo 


afa + 1) OL 
Cele 6 + > (a fB)- 


As a result, (2.20) is negative, which is a contradiction, 
and thus 5 and 8 Ssmould both be positive. 

Now assume , is equal to zero, g then is also equal 
momzero by (2.18). This leads to either or both of the 
following cases. 

a. = © 

This is impossible, since the positive random 
Variable X can not have zero mean. 

b. x =b 

Which contradicts to (2.19) 

Thus, j and g SHOULG. De UWOrhn ser teLly DOSitive as 


required in (2.1). 


24 





Mot CuCsED FORMULA OF PR{X=n } 


As mentioned in the last section, the exact mathematical 
expression of the Pr{X=n} can not be derived by the present 
method of integration. Approximation is then used to esti- 
mate through a closed form of the integral I(a,8,n). 


Hoeou. DY otras Chat 


oo k 
rats = : (-1)* a 
k=0 kK! 
Hence, 
= eg, eaten oes (os 
G1) i@asceam) —=..)— (=1) mr Cor ne 
k=0 0 


Now, consider the Taylor's expansion of 


b -a-n 
Sd army eae 
Se .2) es sa, b i S ays ) Cesap ne eee al) 
; atBt+n KH? k! (atBen)* 
Mii iobyime (3.2) bY Rae ee qos to get 
atB+n-1 ~ k bX F(atntk) 


Zs 





mae, stppose for some values of a and 8, the following 


equation is known to approximately hold 

(ats+n)* = (atBtn)...(a+B+ntk-1) 
je tollows then that 
(3.4) r(atgin)(atetn) = P(atp+ntk) 


and thus (3.3) becomes 


at+g+n-1 mah b“P(atntk) 


(3.5) r(B)b ae kif (atptntk) 


fine into account (3.4), then comparing (3.5) with (3.1), 


it leads to 





oe -l 
Cl + Bye Er(gyp et tne? kota 10,840) 


recalling from the last section, Pr{X=n} is written as 


Pr{X=n } = ea? 8). EL. Geen) 


nib ©" 1(a)r(g) 


Hence, 


b' T(at8)I(atn) 


a a b 
eee Friken} © irrtayPCataen) ‘1 * atgin 


Nc 


26 





HG can be seen that formula (3.6) above is an overestimate 


of Pr{X=n}, i.e. 


= 7 bo’ r(at8)I(atn) le) -a-n 
priten} = air(ayrCarpin) ‘1 * arp? > Fn 


where E. Poo pootur ve quantity, depending on a,B and n. 


mietact, Let 


n k 
(3.7) T! _ pitt) b T'Catntk) = 
° k!T (atBt+n) (atBt+n) 


k — pb rots) br (a+ntk) 


(3-0) T2 = “AIP(Ca)« -RirCatBentk) 


Multiplying both numerator and denominator of Tt by 


Catgtn)...(atBtntk-1) it results in 


(3.9) pe = BUR(ots) béP(atgtn) (atgtn)... (atptntk-1) 
: Pe DT Cae ecree orem 


@eking the difference 6. Gums. 9 )eand (3.0), end noting 


that for 6, = 6, = 0 then for k > e 


bir(atp) b<T(otn+k) , (atptn)...(otptntk-1) _ 44 


Sy, 7 niT(a) k!flatgtntk) (eet 


Meeoucete) | bP (otek) gt 


= i 
n!if(a) k!fCatBtntk) noe (1 + Teun 1} 


afi 





litmus, tor k > 2, Oy increases as long as the numerator 
memrareger than the denominator, then for some k, it converges 
to zero. The speed of the convergence is proportional to 
eer 6. 

It follows that 
k 


(-1) 6, 


m 
i 
ll 1 8 


Hence, for a + 8 large enough, E Powomet ieee o) 1S tnen 2 
good approximation of Pr{X=n}. 
In practice, (3.6) can be applied vastly together with 


some correction. Before doing so, however 6, should be 


k 
expressed more explicitly. 
Consider 
k-1 5 
I (1+ —s—) - 1, 
j=0 o+8+n 
and assume that — is small enough so that 
Claris ce 
k=] 
} MPM ee a, 2 k=l 
420 (1 . ea 1 atetn - Che cas al 130 atR+n 
e (k-1)k 
2(atBt+n) 
hence, 
b' T(a+8)P(atn) (P K(atn). cori: 1) Tena S 


She i: Memeinilat atein) * Gee .(atBtntk-1) 
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and 


b T(a+8)P(atn) f (-1)* b*(atn)...(atntk-1) 


fn ~ 2CatBtn)nil(a)l (aren) — [ou Geuniotn) Ge Otp dtr —) 


(Kav) ks 


Mie factor in brackets does not depend greatly onn. It 


approximately has the common value o for some small n. Thus, 


E : b' T(a+8)P(atn) e 
n 2(atBtn) nilla)rlatBtn) 


now suppose there is an N, such that 


N 
1, Secu) 
n=0 


br (at8)P(atn) 


n!Il(a)ICatB+n) 


b 
(1 + chou al 7 1 


Deer en tenclose to zero for all n > N 


ect, 


br (atB8)P(atn) 


N 
=o - F 
= 2(atBtn nila) lCatB+n) 


Solvineetor 0 in this equation, it results in 


a 





ul b T(atB)P(atn) ia 


o = (S-1) [ 2 acatgen)nil(a)r(aretn) | 


Eweeebeuting © back into (3.10), to get the approximative 


e,» and finally, with the correction, (3.6) is modified as 


bes = b’r(atB)P(atBn) b on 
” 0 } 
2(atBtn) 


The application of the above results will be seen in 
Section V. As is shown, the estimate of Pr{X=n}, using (3.11) 
agrees to two decimal digits, in overall, with the results 


Meme (1.3), in Section I, for b = 6 and a + B = 13.009782. 
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IV. NUMERICAL TECHNIQUE IN APPROXIMATING PR X=n 


In the last section, a closed formula of Pr{X=n} was 
Semived. Although it proved itself a fairly good estimation, 
tne error involved was uncontrollable, and thus the validity 
feepune Droposed distribution of A in Section II cannot be 
mooeed., Hence, numerical techniques of integration are 
suggested to evaluate the Pr{X=n} through the integral 
I(a,8,;n). Several methods have been considered, such as 
Gaussian quadrutures and Romberg's extrapolation to the 
limit [4]. However, the direct application of these methods 
leads to unsatisfactory programming considerations and 
unbounded error estimates. 


Remember that Pr{X=n} can be written as 


n 0 k 
a: _ b TPCatB) k b T(Catntk) 
ens = aay) ae tL) Faeccraeeh rin) 


equivalent to a series of alternate sign terms which converges 
monotonically to zero for large enough k. This behavior 
guarantees that Pr{X=n} can be estimated by summing the 

firses K terms, and the absolute error is no eecey than 


Paewapseluce value of the K + ec Wet. dt. em 


Perna. a ohermen el) 
eich lol ers tat +. 


oh 





It seemed to be a powerful method, however the conver- 
memee 15 rather slow. Moreover, if the estimation is done 
Mmeemmeans Of a computer, the precision will be altered by 
the propagation of the error in repeated subtraction of 
nearly equal terms in the alternating series. 

Thus, a@ special technique iS required which is expected 
to improve the convergence and also eliminate the error 
Peepagation. 

This technique is simply a variant of the Gaussian 


quadruture method. 


A. CONVERGENCE ACCELERATION 


The idea is to approximate the factor aol in the integral 


T(CatB) atn-1 


gate 
(b-)) SN 
pete-lacayr(g) o 


Pr{X=n} = 
1a 


r 


with an orthogonal polynomial pO) of degree s by least- 
squares approximation method [4]. In other words, it is to 
iP aLiaKel pd) so that 
a, 2 
(4.1) E(s) = fs fe - p6A)} di is minimized. 
0 
From the theory of this method pA) can be expressed 


as 


S 
(4.2) pA) = A, She 
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where the {Pj3 is a sequence of ortho-polynomials of degree 


aeeatia Savisfies the properties below: 


Py (a) =]. 
b S. 3; aifke= j. 
(4.3) if Py CA) ip AC aol | J 
0 J OP ole # 4. 
1 De 
(4.4) eee ea PX) ad 
J =F 0 J 


ims mepned also provides a recurrence formula of 
P40) in terms of PO) and P46) which is 


(4.5) P54 6A) = (i - BS) PQ) - cP Oy» 4 


j-l 


where Pe and cs are defined as 


b 
(4.6) Be = = Sf AP.(A)* aa 
J 5 0 J 

and 

oye 

aS Oi) 2 2 

(4.7) C, = el 

@) one j = 0 


For the sake of simplicity, let now denote POA) by P and 
p.(A) by simply Pa: Then, taking into account (4.1) and 


(mee iis) was expressed as 
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mis) = f te -p} ar 
0 S 
b o) S S b s ‘ 
me eect eee d.d. J p.p,dX - 2 2 da. f° -) 
0 j3=0 k=0 ee 0 nige 4=0 - e Psd. 
by (4.3) and (4.4), E(s) is reduced to 
~2b 2: ee 
mes) =el/2 (1 -e °°) = £4,° § 
j=0 J 5) 


Thus, for a value of b, E(s) is decreasing when P. takes 
a higher degree, and as b increases s should be also increased 
in order to keep E(s) as small as desired. 

As an example, some of the values of E(s) were displayed 


imesable V. 
TABLE V 


The least-squares error E(s), for b = 6 


E(O) = «33415549 

CT SL ue ey | 
E(2) =  .20317197 107 
Es) = ee see 092 One 
E(4) = .18656389 107? 
E(5) = .10601508 107" 
Boje -meeeiseco 147 107° 
E(7) =  .15041719 107! 
E(8) = .39915586 1077 
E(9) = °#&.86497233 1Ouem 
BGligy= MMSss653 10°-° 
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B. GENERATION OF THE SEQUENCE {P,3 


eye aetinition, FP, can be written as 


J 


J 
(4.8) P= anc! baie 


where the AJ's are the constant coefficients of the 0 


mas noweto penerate explicitly the {P,}, i.e., to find the 


J 


a's in terms of b. Before doing so, let us express 2a 


B 


5? and ake Ma computable form. 
Taking into account (4.8), (4.3) is then written as 
| b ; 4 
mo Ss. = £ Cf S as Ver EE. 
“ k=0 g=0 0 


J J 3, j.ktgt+1 |b 
= 5 = “ Ms A 
a0) P= wee 0 


J J Area ed Cat 
= r ? Hyp thy |2 
k=0 &£=0 k+Q+1 


Similarly, (4.6) becomes 


i, J 


j SN pgleee 
(4.10) Be = er EF ia an 
j k=0 2=0 0 | 
J J j, j.~kt+£+2 b 
== Er or fk tp? 
j k=O £=0 — ktete 0 
k+9+2 


jee Marwan 0 
eee gD 
j k=O 2=0 ktite2 


Se, 





and finally (4.4) becomes 


4 b _ 
(4.11) d, = = ) A,” f Me aa 
j k=0 0 
3) k 7 i a 
: S, a it a : Cyt bo? + kt} 
= r= 


Paeoarvicular., tor | = 0 applying the results just 


developed to get 


S)=b3 B= b/2; and d= (1-7 )/p 


Then, by (4.5) 


Po al Gon OW A) 


With the expression of P, as above, applying again (4.9), 


1 
Sete wanc (ial) stor 7 = 1, it results in 


=. 7 
Sy [ieee By b/2 and 
d, = -e7°(12/b7 - 6/b°) + 12/b> - 6/bd* 
Hotes liek is generated recursively by (4.5) and 


Shinde 
by means of a Fortran-Formac program [5], which has the 


Sapavllity Of manipulating symbolic mathematical expressions. 


This program also computes S47: Bray and Asay @ivectiy hy 
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fo. 9), 
generated, are displayed 
mone LOllowing is the 


Formac program mentioned 


ALGORITHM 1. 
Step. Read. in 

Read in 

set j = 


beeper 2.) Compute 


Compute 8S 


Compute 
Compute 


SES) Oa 


Cele me and= by 4.11). 


some of the first Pas 


SO 
ijielable VI. 
Algorithm, which the Fortran- 


above is based on. 


(4.5) 
(4.9) 
(4.10) 
(Gi een ale 


Sea ee 


jt PY 
Pye PY 


oa ee 


lisdarane jpy 1! 


GOnREO] SCeD uc. 


With Fortran-Formac programming, 


Ootain Seal 


is known based on (4.5), 


at any desired j. 


given Hoe 


iiwprace1ce., 


e we can 
However, ied > 


the coefficients A's can be 


obtained with Fortran programming alone as follows. 


jJ+l 
=o a) 


ae < m < j-2 


a 
m 

and 
gd J 
AX a6 = Bek. 
j-1 Jj j-1 
7 
A = A 
J j-1 


= oped 72 aie 


= pay —~ cad} 4 J 
jm jae m- 


0 > 


oe 


+ AS 


poet 


Sif 





ABE Vi 
P, and ae HO. i Ses oe 
re jh eRe. 


3 Ee /2 + IWS = b2/20. 


~e7(360/b> + 180/b!' + 30/b2) + 360/b> 


EoD 4 20/b-. 


~e79(16800/b! + 8400/b° + 1680/b> + 140/b") 


+ 16800/b! ~ 8400/b° + 1680/b> - 140/b". 
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But , the AJ's computed with Fortran programming were 
Meenas precise as those derived by Fortran-Formac, since it 
was dealing with a large amount of computation, while the 
AJ's are evaluated directly from symbolic expressions in 
Fortran Formac programming. However, the differences were 
negligable when j was less than 10. 

iews, taking account of (4.8) and with the AJ's found as 


above, P. is now written as 
S J 
ea 12) p = 1 Scape 


Now, let P(n,s) be the estimate of Pr{X=n}, approximate 


e with p.. Taking account of (4.12), 


S j : b 
P(n,s) = += Ed, f Ae ee) aac) 
n!b PCa) PCB) j=0 “k=0 0 
‘tees j 
Further, let ty ee AY — (hens 
el Ss 5 : 

Beeman ee ieee creecc Ertl) 

CR SS aera — weg Kk > FCar84ntk) 


and P(n,s) above was computed by means of an iterative 


method, based on the following algorithm. 


ALGORITHM 2. 


Seepea, cet j = 0 


n 
_b T(atB) ,o PCatntk) 
Let P(n,s) = n!if(a) a T(atB+ntk) 


Se] 





DPiepee. oC.) = J] + 1 


br (atB) J k T(Catntk) 


E j 
Compute TCJ) = rFtay cig *k © FCatBmtk) - 


coo Penis ea b(n es) + TC) 


CCE ei) | = 2 Go to Step 2 


ive iec tr ees LoOp 


It is noted here, that Step 2 was called the J + ve 


iGeration. The stopping criteria in Step 3 is only for the 
purpose of fixing the idea, it was replaced later by a more 
appropriate one discussed in the next subsection, when the 


error in P(n,s) is discussed. 


Cee kROR ESTIMATION 
Let e(n,s) be the error committed in approximating P(n,s) 


fomrrixX=n}. Thus, 


as) ExGamice) Potten) —wbm.s) 


b 
E Lee ee 4 AMPA Eben) EN (e™ pg) ad 
née Bi 


Now by the mean value theorem of integration, and by 


atn—-1 y 8-1 


foringe that 2 (n=) does not change sign in [0,b], 


Giewetierc exists a A*eclOpb], which depends on Mee (bak 


and on ep (A), SOmeEmMat 
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= T(atB) 


Wy Baeave cate 


b 
(4.14) e(nys) = fe7*"-p_(a*) oe ee yak 


In general, such a A* in (4.14) is not known, hence e(n,s) 
cannot be found. However, by using the Cauchy-Schwarz 


inequality the bound u(n,s) of the error is derived as follows. 


ail 
b + 
feim,s)} < u(n,s) = - ee ae y2Catn-1) py) 28-1) 42 
Ib PCa Cs) 0 
i 
b + 
c s (es 1 BONE 
0 S 
hecalling that, 
b S 
Hos) = ff ta 6 ye dA = : “le Gs r ac Sa. 
0 S j=0 J 9 


and that 


4 9) 
2(atn-1) 2(g-1) _ , 2(atgtn-1)-1 [(2a+2n-1)F (28-1) 
fo Co-d)- dA = b T(2at opton-2 


0 
then 
N-Br(ats) (1 )P(2p-1) 3 
ome cl Cons cL 2aten-1)T(2p—-l 2 
4-15) uns) = orevayrce) F(Zateeren-1) «= SS?! 


On the other hand, the error e(n,s) in (4.13) can be 


alternatively written as: 


4y 





b 
om, S ) a ae if yoemat ny) Pte an 
n!b ccnp 0 
b 
T( ats) | atn=-1 B-1] 
_ r (b-A) Ol AG 
~ i 
nibo ’ “P(a)T(B) 0 : 


Now, applying the mean value theorem of integration to the 


two above integrals, and let 


We =_ Peete) Pi yetn=1 (44) 8-1 ax 
Gee, 
nib” T(a)T(B) 0 


_ be? (at+8)T(atn) 


nil(a)lCatb+n s oes 


r 
- =e 


where A, and dy SGOn bl 


Q 
-i 
Modine and subtracting to (from) (4.16) the quantity e 1 Vo 
we have, 
a “4G fia! 
C4. 17) e(n,s) = fe _ p,(A,)} V + fe -e } V 
=} =) =) 
=e) t= P(n,s) + fe 6 a V 


Further, assume that do ~1> then (4.17) becomes 


a “AS 
Cae) e(n,s) =e V - P(n,s) 


4a 





It can be seen that (4.18) is a good estimate for e(n,s) 


only when the above assumption holds. 


The A, in (4.18) can be found by solving 


al 


en 19) 


as, 


GAY) = p,(A) as Ein 5). = 0) 


Based on the above analysis, algorithm 2 is now modified 


step l. 


SUED Cc. 


step 3. 


Step 4. 


However, 


set j = 0 


br (ats) o I'(Catntk) 


Set Pin,s) = Tatta) *o Tratéentk 


Dees ject a 


= br atB) : tJ pk TCatntk) 
mea) bet T (atgtntk) 


Compute T(j) 
POC 5S) Ge) as Bees) 
Compuce Ulm. s) 

Wie SOS)! te Stop 


iE AU Gol, Si) mene Gente otep ec 


fore the purpose of comparison, the approximate 


moneeoeas an (4.18) is also computed, and it is suggested 


that (4.19) be solved by Newton iterative method which is 


Hecate or Gnas problem {4]. In fact, f(A) is continuously 


differentiable in [0,b] as the Newton method required. In 


addition, the derivative of f(A) is easily computable. 
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taking into account (4.12), then 


S j : 
Pe) = op eel eee Eins) 
j=0 k=0 
i lollows that, 
S J 
(4.20) = FOt) = es ae . ,k-l 
j=l k=l 


Moreover, (4.20) evaluated at A 
r 


1 does not vanish as Pp 6A) 


approaches e- 


60) 


On the other hand, a starting point for the Newton 


meeravaons is nicely provided by: 


(0) 
na 


_ Eanes) 


(4.21) 


In such conditions, the Newton methods are very efficient, 
i.e., converge rather fast (quadratically) to i: 
For a detailed treatment of the analysis, the reader is 


referred to Appendix A. 
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V. AN ILLUSTRATION 


Consider the following example. 

A sample of 621 white male children is observed, when 
they were from 8 to 11 years old. The data reported here 
mma part Of a large study of childhood accidents conducted 
by the State of California Department of Public Health. 
Subjects to study were those children whose families sub- 
scribed to the Kaiser Foundation Health Plan [6]. 

The number of injuries and the observed frequency were 


tabulated below. 


Observed Frequency 


X 

0 240 

il 192 

2 107 
pie 


| 
and more k 
For purposes of comparison, let us first apply the 
Negative Binomial to these data, i.e. assuming that X is a 


Poisson variate with parameter A, and A itself is Gamma 


distributed. 
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A. THE NEGATIVE BINOMIAL MODEL 
iw swe DUG ToOnson i 
Using the standard method to estimate the mean and 


the variance of X results in 


decuins 


> | 
i 


and 


See wl. oc 


Memee, by (1.4) and (1.5), 


m= 1.13 


and 


2 
I} 


bee mitic wk /m = 2507 


iio bowcutaan frome (!.2). 2A) is Gamma (3.24,2.87). 


The mean and the mode are respectively, 


ELA] ais des 


d = .785 


The curve representative (C) of g(A) was shown in 


Heapte Vil. 
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ae emunmcondiutonas Probability of K =n 
With the values of k and as above, and by (1.3), 


Pr{X=n} is 


3.24 +m) is se ee 


meen, for n = 0,1,2,3,4,5,6 and more: 


n Pr{X=n } Expected Frequency 
0 . 38063 236.37 
if os > 197.38 
2 ~17374 NO, Giele 
3 .078 36 48.59 
4 .O3149 1S Pas: 
» Bon wilyas (far sie, 
> 6 .00627 3.89 


x" = 1.179 with 4 degrees of freedom 


152 < p < 90% 


Seer lCA EON MOn SECIUTONS Til, Til, IV. 
Paleo ripurion of A 
Assume now g(dA) is of the form (2.1), and assume 
further that the probability of X is greater than 6 is small, 


then by (2.9) 


ail 





| 


Se NI eI AI |, ( 
Cae oa a 


(J) SAGean tg 
nen: 
40i 





een 
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Merollows by (2.17) and by (2.18) that 


a = 2.440642 
gB = 10.569140 
and. 


a + 8 = 13.009782 
hence, 
SO cies pee 26a: 


Ome sls 10°” (ON 


with mean and mode respectively, 


1 


ALA) = 1.13 


RO Sais 


The curve representative (B) is shown in Table VII. It 
is noted that (B) and (C) have the same mode, Seen O15 
Pa UoemUnCOndt Lona —Fropabi lity 
a. Using Closed Formula. 


Substituting the values of a and 8 found above in 


(SM, 


8 
cee oo oa To .iooke + pn) 
cate = n! P(13.009782 + n 


6 y-2.440642—n 


(1 + —————_—_—_ 
oe e0o fee tn 


1 





Then, for n = 0, 1, 2, 3, 4, 5 and 6, the rough estimates 
Of the probabilities of the unconditional of X equal ton 


are the following 


n Pr{X=n} 
0 39629 
i 5 SCL 
2 .18628 
5 .18628 
4 .08685 
5 .01338 
6 00463 


Now applying the correction to the above results. 


memes. 12) and (3.12), 


So = 1. 053.3 
S-l = nO 33 

and 
a = F901. 


Thus, the results should be corrected as 
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n Corrected 


0 . 38091 
iL .31049 
2 ee ere 
3 me Op a 

SUSIE S 
> eC) ace lea 
6 00414 


teecameme seen tLnav With the correction, i.e., 
using formula (3.13), the estimates of the probabilities 
of X equal to n, can be compared within two decimal digits 
wemone results found by the Negative Binomial model. However, 
the significance of these figures is uncontrollable, hence 
there is no use to perform a goodness of fit test here. 

b. Using Numerical Technique 
Hii iiveacdvance r= 2156 nee fHeneatcer Ll 


iterations P(n,s).converge to the following figures. 


n En. Ss) Ge 4) Expected Frequency 
0 . 383299 2.0 10! 238.03 
1 . 313539 D0) ial 194.71 
2 173662 Ait go! 107.84 
3 .0795 39 142 ao! 49.39 
y .032189 $.2 oe 19.39 
5 .011869 2.8 107° ies i 
6 OHNO ih ge 51) 


aul 





Based on the above results, and by noting that 
fre sum Of P(n,s) from n = 0,1,... equals to one, a chi- 
square goodness of fit test is performed, the results are 
y° = 1.035 with 3 degrees of freedom and 75% < p < 90%. 

foicsemnercecdmoaac the values of the P(n,s) and 
Mem. s) are rounded off to the significant digits. For the 
convergence and the actual values of P(n,s), u(n,s) and 
e(n,s) the reader is referred to Appendix B. 


The summary of the results obtained in this 


section are shown in Table VII. 


By 





On WW FP WW PO KF © 


Note: 
(i) 
Calat)) 
Casa) 
(iv) 
(v) 
(vi) 
(vil) 


el) 

. 38063 
e735 
» 
.07826 
-03149 
mo 7,.6 
.00627 


(ii) 

. 38091 
.31409 
p22 
SIO Gal 
- 03279 
see Ul 


00417 


TABLE ViliI 


Gait) 

. 3832990 
5 SUES S15 
.1736619 
-0795389 
.0321885 
.0118698 
.0059066 


SUMMARY OF RESULTS 


(iv) 
240 
192 
Leys 
ae 
iby 
9 
4 


Pr{ X=n} (negative binomial) 


(v) 
236. 
dS ae 
LON: 

48. 

19. 


(vi) 


2362 
04 


we 


108. 
90. 
Zs0le 


Pr{X=n} (closed formula 2 with correction) 


Pr{ X=n} (numerical 


Observed frequency. 


Expected frequency 
Expected frequency 


Expected frequency 


BS: 


technique) 


(negative binomial) 


(numerical technique) 


81 
M2 
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ae 
ol 


(vii) 


238. 


194 


ve 


ial 
iO ie 
Hg. 
19. 
37 
. 66 


39 
98 


(closed formula 2 with correction) 





Vale CONCLUSIONS 


As is shown in Section V, an application to the real 
situation of the new model is presented. In it we assume 
Maas the data 1s Poisson distributed with parameter A and 
the parameter itself is also a random variable with density 
iumeri1on of the form as specified in (2.1). The validity 
of the model is then made by performing a chi-squares good- 
Meso of fit test. As the results, ve = 1.035 with 3 degrees 
of freedom and 75% < p < 90%. 

On the other hand, modelling the set of data as Negative 
binomial model, and using the same validation procedure the 
results are sf = 1.179 with 4 degrees of freedom and again 
752 <p < 902. 

Comparing the two results, and referring particularly 
to the expected frequencies, it can be seen that those of 
the new model are closer to the real data. However, because 
the distribution specified in (2) Shes 3 parameters, i.é., 
i paramever more than the Gamma, causing 1 degree of freedom 
less in the chi-squares goodness of fit test the results do 
Mewemciiter sre@nificantly. Availability of an extra parameter 
is one of the major advantages of the new model in regards 
to the modelling of any real situation. However, using the 
maximum of the observed values of X in a single period as 
an estimate of the parameter b in Section II, is not an ade- 


quate method. It leads to difficulties in extending the 
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mecel in order to predict the distribution of X in the 
next period. 

Moreover, the exact mathematical expression of Pr{X=n} 
@emney be derived, hence the Knowledge about the statistical 
meeperties of X is limited to the mean and the variance. 

iimeoection Ill, & closed formula of Pri{X=n} is then 
Gerived. However, it can be applied only in cases where 
the sum of a and 8 is large enough. 

Finally, a numerical method for estimating PriX=n} is 
presented in Section IV. It can be seen that this is a 
powerful method when high accuracy is desired. However, 
it is also complicated and can be best used by means of a 
high speed digital computer. 

As recommendations for further investigations, the 
following are considered as potential topics: 

1. Determining the exact expression of Pr{X=n}. It 
is believed that such a formula can be obtained by advanced 
mathematical analysis. 

Ce Hetemmdnine an adequate estimation procedure for 
the parameters involved in the prior distribution of A as 
Spceowmued sn (2.1). 

SECON CInMIngG wath aAnewexmtensdon Of this model to a 


bivariate or more general, k-variate model. 
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APPENDIX A 


COMPUTERIZING THE CALCULATION OF P(n,s) AND 


DETAILED TREATMENT OF THE ERROR ESTIMATION 


The following shows the development of the modified 
algorithm 2 in Section IV into details to computerize the 
calculation of the P(n,s) and the error involved u(n,s) as 
well as the estimated error e(n,s). 


S Be ane 


Ome 0) 0) de 
1? By were already input. We now wish to compute P(n,s), 


Suppose the values of a,8,b and Pos qd 


d S 


1? 
Menamcand the corresponding e«(n,s) for n = 0,1,...,N. 
since the computation is done by A on eee method, the 
Swepoing criteria should be determined first. 

Let r* be the maximum of the u(n,s) for n = 0.1...,N 


th 


at the s iteration. The iteration is then stopped when 


r* < r, where r is a small quantity determined in advance. 


Step 1. Initialize all the P(n,s) to zero. 
-2b 
) 


Meso evitere sls tyne Subscript of the Bac 


Seu kisi) — .5(1—-e 


UCC mEU Omnium. 55..65 and note that s is the actual 
Valve enn) ne this step, then 
Compute Pas Ss Bes d. DemGes C4. Oso (4. 10). 


and (4.11) 


J 
ae 


Set E(s) = E(s) - (ad 


Compute t 


o 
as 
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meets eee Or tf = 0a... N and for j = 0,1,... 


b’r(atB) k T(atntk) 


jo 
j 
nif(a) “5 “k > PCatptntk) 


Compute T(j) = 


wee 1) SotmuGaas) samen.s) + T(3) 
Compe uin.s) by (4.15). 
Cempuirewewm.s) by (4.18) for 37 > 0 
tersceby Solving for hy Lie: 


(a.2) 9 £0) = p(n) - SRS? = 0 
As mentioned previously (a.2) is solved using Newton 
iterative method [4]. The idea is to generate the sequence 


,62)} where i = 0,1,... by the iterative formula: 


(m) _ £(aS™)) 


, (m+1) 2 
£1(, 6M) ) 


(a. 3) X 


starting from , 60) as defined in (4.21). 


Furthermore, if the sequence so generated converges to 


some point €, i.e., for some inter L, we have 


(a. 4) , (L+1) - 6h) ie 


so that (a.3) now becomes 


B= 8 - tS 


out 





then it is immediate that f(&) = 0, and that & is the 
desired i: 

Peotomanary, by (a.3), £(A) should be continuously 
differentiable and its derivative with respect to A should 
mew vanish at € = i: However, referring to the analysis 
done in Section IV, it can be seen that f(£) satisfy both 
of the above conditions. 

On the other hand, starting with 60) ao Orman 4.22) 


the sequence eed will converge to A. quadratically. The 


i 
Peeot is as follows. 


Consider, 
(a.5) h(A) = dA - 


Taking the derivative of (a.5) with respect to id, 


Blas 
(p_(A) - ——) f"(A) 
(2.6) i ee 
Ee eny | 


ewe y Noting that, in the vicinity of A (576) as 


le? 


small with respect to 1, and it vanishes at A then making 


1? 
use of the discussion in page 57 of the reference, it follows 
that the iteration converges quadratically to ds: 

But, the integer L as in (a.4) is not known, and usually 
is estimated by L* in advance. Hence, x, may be found 


exactly or with some tolerance within L* Newton iterations, 


ef 160) is close enough to i: 
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Based on the above argument, we may then stop the 


searching process for A, if for some m the following conditions 


1 
are met. 
(a.7) j Om) = ym) < $ 
(a.8) ie(aS™)| < @ 


Grelesnovedstm@at m < L* and that ¢, @ are small positive 
quantities. 
The above conditions are checked by associating them 


with the binary function Flag defined as below 
Ot 62-7.) or (a.30) are satisfied 
Flag = 
1 otherwise 


Then during the course of searching for a i af eed 


lew 
comes out with the value 1, either one of the following 
corrections have to be made. 


1. Increasing L* 


(0), (LD 


Ae oUDSGlvubaae A by A 
The corrections may be repeated until Flag = 0. 
Summarizing the above analysis in the algorithmic form 


we have: 


a2 





TH Ouee& |G 0 esee ees 

Compute , 60) ye Coe 1) 

Solve for hy Paeicase ) 

Find the value of Flag, make the correction 
it Hees = | unter! Klap = 0, if 60) is close 
enough to dy: 


Compute e(n,s) 


Step 4. Check if u(n,s) satisfied the stopping criteria 
Was oh Oe e si. sy Ns 
hiGees 2 Ss UOp 


hmMewciaee re Updave j by j + 1 then go to Step 2. 


The application of the above results is seen with the 


Mmeiiseration in Section V. 


60 





APPENDIX B 
TABLES 


mae COMvelivse et thts appendix are for the purpose of 
miustrating the convergence of P(n,s), for n = 0,1,...,6 
as in the example presented in Section V. This is also 
known as a direct application of the results developed in 
section IV and in Appendix A. 


iiemwvoatues ot P(n,s), u(n,s), e(n;,s),;. A f(r) and Flag 


1S 
at the j th iteration (j=2,3,...) are tabulated in tables 
Mmeonoered from 1X to XVIIL. 
The following abbreviations are used. 
P(X=n): 124 Goins) ies aime Cel eee oe 
BOUND: ie boundect the error, 1.e., Uuin.s). 


EST. ERROR: The estimated error, i.e., e(n,s). 


oC Oe The value of dy: 
FXO: f£(d,). 
JDM iy Mee Flag. 


It is noted that the accuracy desired is achieved at ll 


2 


iterations, with r = .156 10 ~, the tolerance error in dy 


~15) 


and in f(r.) less than Ue? (d = 6 = 10 and the maximum 
number of Newton iterations allowed (L*) at 15. 
The generated tables shown in this appendix and all compu- 


ter work presented in this thesis were done utilizing the 
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IBM 360/67 computer at the Computer Facility of The Naval 
Postgraduate School, and programmed in Fortran-Formac and 


Mmierortran G using double precision. 
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PRtx=-0),@CUNC 
feole ERROR 
XC,FXCyIFLAG 
FR(X=1),ECUNC 
EST. ERRCK 

XC sFXC si FLAC 
PR(X=2),@CUNC 
EST. ERRCR 
XC,FXC,IFLAG 
PR(X=3),eETUNC 
EST. ERRCR 
XC»FXOCyIFLAC 
PR(X=4),BCUNC 
EST. ERRCR 

XC ,FXO,IFLAC 
PR(X=5),ECUNC 
Gore. ERRCR 

XC ,FXC ,LIFLAC 
PR(X=6),ECUNC 
EST. ERRCR 

XC ,»,FXC,LTFLAC 


R > ic 
GREATER THAN 


ACCLRACY CESTREC HAS NOT BEEN ACHIEVEC, 


R 


TABLE IX 


ITERATION # 2 


O22 252540 


0.11256032D 
0 .37899285D 


0-14/735313D 
0.25150315D 


O ei 7509920 
0.13423654D 


0.20389942D 


0 .615922828D0- 


0.22/718605D 


O25 505 9300- 


0.24 /88669D 


Cmep orc lLol— 


0 .26640943D 


ere 


Cl 
CO 


Ol 
CC 


Ol 
CC 


Ol 
Cl 


Ol 
Ol 


01 
Oz 


Ol 


Ores. 1 3160) 700 
= Cero go e530-0 | 
0.138/7/7880-16 
O.2451741900 OCC 
~OstelLlLoocgD: C0 
0.0 
CellE47110D0 OC 
~C.L1L07€249D0 OC 
C.1SE77178ED-16 
UslLUsoos21D CC 
=Gel026/2410-01 
Gece fi (830—=16 
Cee 2i olor pol 
=Us 201290 14D-01 
Oetocii(ec0-1 6 
Os240c(So7cu-O1 
aUelaocO4C90-Cl 
O, LOC ermse cD=1¢€ 
C.10024512D-01 
=O 36256 1400-0 2 
Gee 7 4sG2SIR=16 


0-2451719C30 OC 
= O0.15600CO000D-05 


ITERATICN CCATINUES.. 
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PR(X=0),8CUNC 
EST. ERRGR 
XCsFXC,IFLAC 
PR(X=1),ECUNC 
EST. ERKRCR 
XCsFXC,IFLAC 
PROX=2),PCUNC 
EST. ERRCR 

XC ,FXC,ITFLACG 
PR(X=3),@CUNC 
EST. ERRCR 
XC»FXO,TFLAC 
PR(X=4),8CUNC 
EST. ERRCR 
XCsFXC,IFLAC 
FR(X=5),ECUNC 
EST. ERROR 
XC,FXC,IFLAG 
PR(X=€),ECUNC 
eels ERRGR 
XCsFXC,IFLAG 


R* ies 
GREATER THAN 


ACCURACY CESIREC HAS NOT BEEN ACHIEVED, 


R 


TABLE X 


ITERATION # S 


0-41672315D OC 


0.10631004D Ql 
0.3€144446D C00 


Orls9399430 Cl 
0-21239926D QO 


0.-16813740D Ol 
0.98482140D-Ci 


0.19330664D Cl 
Oee9169598D-01 


ez t 205/90. G) 
0.13 7758050-Cil 


YecsleloUzZ DOL 
0 .43490496D-C2 


06252795640 Ol 


0.9902C114D-01 


SUeile e3e240—01 


OL 40 1TECCGZ2D-—1é 
Oetwe45solbD OC 


=O2082C21LZ00-C1 


Os 25CSCGGD— 16 
Octo dor C54an-—Cl 


-0.5805C47/D-O1 


Ce 2040695 10-1 
0.4455S8/77D-01 


=—Qezial(ssslil>Uur 
=O. L2twereepa1 é 


0.22/712345D-01 


HSCerOrces7¢oD=01 


O.13877788D-16 
Gel 3655730-01 


-€.29€0C 2500-02 


O227 tee otoU-—LC 
Cet al (SsiSp=0< 


=CeoUG Si o050—-05 


0.0 


0.1054S960SD 00 
= 0.15600CCOQO0CD-05 


ITERATICN CCATINUES.. 
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PR(X=C),ECUNC 
ols ERRCR 
XC,FXC,ITFLAC 
PR(X=1),BCUNC 
EST. ERFRCR 
XCy,FXOC,IFLAG 
PR(X=2),BCUNC 
EST. ERRCR 

XC eFXO,IFLAG 
PR(X=3),8CUNC 
EST. ERRCR 
XCsFXC,LFLAC 
PR(X=4),2CUNC 
EST. ERRCR 
XC,FXC,IFLAC 
PR(X=5),ECUNC 
EST. ERRCR 

XC yFXCsTFLACG 
PR(X=€),ECUNC 
EST. ERRCR 
XC,FXC,IFLAC 


R ¥ hs 
GREATER TRAN 


ECCURACY CESIREC HAS NOT BEEN ACHIEVED, 


R 


TABLE XI 


ITERATICN # 4 


OeosgeczcslSpD CC 


0.10034698D Cl 
0.32485 /736D QO 


Q0.131/70308D Cl 
0.176944670 OC 


0.15893C056D Cl 
Oni SoicetocD—C 1 


O.18272949D Cl 
0 .30485588D-O01 


O.203636¢01D OQ] 
O-lL07@25/7CS0-Cl 


O.22206769D Ql 
0. 44992633D-02 


0.23834%402D Cl 


Oreo (66 5S20-01 


~Osecel(ooz0-0] 


Cre 025760C10= LE 
Oe 35914645D-0 1 


wOoecoslde0'0-O1 


C.SO72E1LI4D-17 
Goa eveles Given bi Opt 


Testi lL OlcColl-O0 2 


Ue2htosetoU-1é 
0.15164€13D0-01 
We5991032710-03 
0.97144515D-16 
0.77447542D-02 
0.21573024D-02 
Gee2zcec255(D=l¢é 
eos o25ap— 02 
OelodaTessu— 0 2 
Oe cece noo 1é 
C2147238730-02 
Ue 5) -62/6D-03 
OelS9452200-1¢ 


0.359 7464520-O01 
= O0.-156CCCOCCD-05 


ITERATICN CCATINUES.. 
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PeGxX=Go,ECULAC 
EST. ERFRCR 
XCsFXC,ITFLACG 
PR(X=1L),ECUAC 
EST. ERRCR 
XCyFXC,IFLAG 
PR(X=2),ECUNC 
EST. ERRCR 
XC,FXO,IFLACG 
PFR(X=3),2ECUNC 
EST. ERROR 

XC sFXCs,IFLAC 
PR(X=4),28CUNC 
ES. FREER 
XC,FXC,LFLAC 
PR(X=5),8CUNC 
EST. ERRCR 
XCsFXC ,IFLACG 
PR(X=6€),8CLNC 
EST. EFRCR 
XCsFXO,IFLACG 


R* es 
GREATER TRAN 


1 


ACCURACY CESIKEC FAS NOT BEEN ACHIEVED, 


R 


TABLE XII 


ITERATION @ 2 


0 .38539646D CC 


0.968632340 CO 
Q0.31364717D QO 


QO.127831480 Ql 
0.17234396D CO 


0.15520382D Cl 
Oath ooo 7 rob-—E 1 


Ocelvses3650 Ol 
Ose (CCOSD-Cl 


0.-201837060 Cl 
O.-116401170-Cl 


Oe222124790 Cl 


0.3559688130-C2 


0.24092410D Ol 


Oneaese2i (0-02 


=O. (546e0>530-02 
—-C.520417040-1? 


OelOrgSss7 10 er 


= Cobo sec sone G2 


OPW SB bor) Oba 
O-f27E711560-02 
Ceo So PeLl24)=0 2 
0.433€8C87D-1E 
O0.42€SSEC3D-02 
Oesere 205 60-02 
O.3 az 0230-16 
O.214€4258D-02 
O.1665€0150-02 
C.5945914840D- 16 
aoe Cr logD-Cs 
Gsece220C30-03 
Orr cec2oe7U-— LC 
C.4137€932<D0-02 
Gee 22624820505 
Oszs4SC2590-—l¢ 


QO.101095714D-01 
= O.-L56CCCOO0CD-O5 


ITERATICN CCATINUES.. 
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PR(X=C),ECUNC 
EST. ERFRCR 
XCysFXC,IFLAC 
PR(X=1),2CUAC 
EST. ERRCR 
XCysFXC,TFLAC 
PR(X=2),B8CUNC 
EST. EFROR 

XC eFXC,IFLAG 
FR(X=3),28CUNC 
Soe CRRER 
XC,FXC,IFLAG 
FR(X=4),8CUNC 
Ess ERRCR 
XC,FXC,IFLACG 
PR(X=5)eECLAC 
EST. ERRCR 
XCysFXO,IFLAG 
PR(X=6),ECUNC 
EST. ERFRCR 
XCysFXC,IFLAC 


R* Is 
GREATEF TRAN 


ACCURACY CESIREC HAS ACT BEEN ACHIEVEC, 


R 


TABLE XIII 


ITERATION # 6 


0.38342695D 00 


C.S$58605490 CC 
0.31224038D OC 


0.12740148D Q1 
0.17235804D CO 


Grab 5109650 101 
0./9420€33D0-G1 


QO-1813511i6D Cl 
Gme21835230=C1 


0.204653790 Ol 
O02 1138S2659D-Cl1 


0.225948410 Ol 


0.407474920-02 


0 -24540726D Cl 


One 2e7 2 cOol-= 02 
Ga2zo LET S50D—06 
e893 4S320D-16 
C.240SS208D-02 
Osho 0 2520-0 2 


sees 2G2451D=1¢e 


One t 2467 550—C0¢ 
0.14157378D-02 
G2525¢2E€CS0-—16 
0.1017&784D-G2 
0.618S648/2D0-03 
Ge4Ge245540-1€E 
Ceres 650-03 
0.1290/7845D-03 
C.86/73E1/4D-18 
On2S6e71620-03 


—-G.26/717150D-C4 
=O.9 235591 voU =O 


0.98€24278D—-04 


—Oes726e 263D=04 


Cat 4o5126220-15 


Oee4eSsecelp-—O2 
= 0.15600C000D-~05 


ITERATICN CCATINUES.. 
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PR(X=C),ECUNC 
EST. ERRCR 
XCsFXCeIFLAG 
PR(X=1L),PCUNC 
Bole ERRCR 

XC yFXC,ITFLAC 
PR(X=2),BCUNC 
EST. ERRCR 

XC sFXC,ITFLACG 
PROX=3),BCUNC 
EST. EFRCR 

XC sFXCs,IFLACG 
PROX=4),ECUNC 
EST. ERRCR 

XC eFXC giFLACG 
FR(X=5),ECUNL 
EST. ERRGR 

XC sFXC eI FLAC 
PR(X=€),ECUNC 
Ets ERRCR 
XCsFXCsI FLAG 


R* ges 
CREATER THAN 


ACCURACY CESIREE HAS NOT BEEN ACHIEVED, 


R 


TABLE XIV 


ITERATION # q 


O.S63299712D 


0.95821145D 
0.31249161D 


O.127/70189D 
0-173648390D 


0.15628372D 


0.7955C6890- 


0.18200144D 


0.32201490D- 


0.20519853D 


O.11876448D- 


QO.22616411D 


0.40594618D- 


0 «24515309D 


CO 


00 
CC 


Cl 
CO 


Cl 
Cl 


OL 
Cl 


Cl 
Ol 


Ol 
02 


Ol 


Ono ts 46 o20 202 
0.28859540D-02 
-0.-384891/770-1? 
O.492C71S 10-02 
0.40242€8 30-03 
Oss eee cecUs IG 
O.358 504770203 
Ce 249 2269003 
C.25424541D-16 
0.2102 /C66D-03 
=Cetoioec lau O04 
C.10842022D-1 ¢ 
0.101 226%6D-—023 
-0.6442€S50D-04 
0.241234958D-16 
C.489555S520-04 
#=6.3601(2277D-—04% 
C.99€ZE015D-16 
Us 20ee5 3230204 
“Q.-137C91410-04 
Ose> es se6510-15 


0.498C/1909D-03 
= 0. L56COCOCCD-O5 


ITERATICN COATINUES.. 
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PR(X=0),28CUNC 
mots ERRCR 
XC,FXCyIFLAC 
PR(X=1),ECUNC 
eel. ERRCR 
XCsFXC,IFLAC 
PR(X=2),8CUNC 
mole ERRCR 
XC,FXC,IFLAG 
PR(X=3),8CUNC 
Bots ERRCR 
XC,FXC,TFLAG 
PR(X=4),8CUNC 
EST. ERRCR 
XCyFXC,IFLAC 
PR(X=5),6CUNC 
Esme. ERRCR 

XC sFXC,IFLAC 
PR(X=€),cCUNC 
bon. EFRCER 
XCeFXC,IFLAC 


R* iS 
GREATER TEAN 


FCCURACY CESIREC HAS NOT BEEN ACHIEVEC, 


R 


TABLE XV 


ITERATICN # & 


Geo eec7e4 7D 


Oaoocio13 1D 
Jono > >> 60D 


Oeiz 7607190 
0.1736641 8D 


0.15636408D 


0.195414800- 


0.18199294D 


One a rso50Z0- 


0.20507474D 


CC 


OC 
CO 


Cl 
CC 


Ol 
Cl 


Cl 
Cl 


Ol 


OL1L3869791D~Cl 


Oiec2 99971250 


Cl 


0.402 WG559D-C 2 


0.244€6248D Ol 


Uneo2e 8 1560=04 
0.7248109/70-04 
C.45424E4/0-1€ 
C.90775300D-04 
0.255712920-04 
C.954CS 79 1D-1LE 
0.65338715D-04 
sUs 042 2c560—04 
-€.46112414D-1LEé 
0.3834C 7ESD-04 
=Ueco8sCisaD—C4 
Caoloie 2570-15 
0.19542441D-04 
—lekcecs2Oo2C)~04 
Oesedeze2es 10-1 5 
Omae2e22210—05 
=Uecoe cot 130-05 
Unameserso 7 0—15 
0.371525440-05 
0.43041158D-0€ 
CoCo cerscsU-1li 


0.90 715299 1D=04 
= QO.L560CCOCOD-05 


ITERATICN CCATINULES.. 
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fren —Cd,ECUNC 
fei. EFRCR 
XCsFXC,IFLACG 
PR(X=1),2CUNC 
EST. EFRCR 
XC,FXC,IFLACG 
PR(X=2),ECUNC 
EST. ERRCR 

XC sFXC gIFLAC 
FR(X=2),8CUNC 
EST. ERRCR 
XCyFXC,IFLAC 
PR(X=4),ECUNC 
EST. ERRCR 

XC »sFXC,ITFLACG 
PR(X=5),ECUNC 
Eis ERRER 
XCysFXC,IFLAC 
PR(X=€),ECUNC 
EST. EFRCR 
XCeFXCyIFLAC 


R * Es 
CREATER TRAN R 


FCCURACY CESIREC HAS NOT BEEN ACHIEVED, 


TABLE XVI 


ITERATION # 9 


Uesescso8lLb C0 


C.95892C26D CC 
0.31353855D CC 


0.12781 740D Cl 
Q0.173662270 CQ 


Ol C3o257) G1 
Osi J2290L10-Cl 


QOe-l181961820D Cl 
0.32158494D-Ci 


0.20503695D Ol 
GetVees(340-Cl 


Omacoo90925)D Cl 
0-405/722990-02 


0.24485899D Ol 


Omrooee ste 04 
C./72014410-05 
Oncol ce ciel S0- 1 / 
0.147&?73460-04 
~C.688&44/7E8D-05 
GO. 120 Ces 2p] LG 
Q0.10€427120-04 
—=Ose 04 sa oo1L0—05 
CelEveoOs340-—1e 
0.6245/73210-05 
Oe corce= 5 10—us 
Uae ores ocp- TC 
0.31824 7420-05 
C.535171488D-06€ 
—~CeoravecurO— 1 7 
Caras sa 5130 —05 
0.855€7C1L0D-0€ 
Gs foe lazcop-1eé 
Oeo0 522241 0-0C6 
0.4384S225D-CE 
= UsloesseCGD— 15 


0.147873455D-04 
= Q.15600COCCD-05 


ITTERATICN CCATINLES.. 
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FR(X=0),eCUNC 
EST. EFRCK 
XC,FXC,IFLAC 
PR(X=1),8CUNC 
Bote EKRCK 

XC ,FXC,IFLAC 
FR(X=2),28CUNC 
Eolee ERRCR 

XC ,FXC,IFLAG 
PR(X=2),ECUNC 
Beds EmRCR 
XC,FXC,IFLAG 
PR(X=4),ECUAC 
EST. ERRCR 

XC »FXC,IFLACG 
PR(X=5),ECUNC 
eowme EFRCR 
XCyFXCsITFLAC 
PR(X=6),ECUNC 
EST. EFRCR 
XC,FXC,IFLAC 


F > 0s 
GREATER THAN 


R 


TAB EE XVII 


ITERATION # 10 


0.383299C6D CO 0.204332300-05 
~C.15S1C&69D-0€ 

0.95894C18D CO -C.242£75610-1€ 
0 .31353853D CO 0.217€81C040-05 
-C.L7G9EL116D0-C5 

Orliz7gls7sD C1 0.11Z2€S229D-15 
0173661980 00 0.15€€8&358D-05 
-0.505562540-CE 

0.15634810D Ol -C.28111541D-1€ 
0.795389190-01 C.91541¢653D-Cé 
0.400435S0D-0€ 

O-NsIS538i40 01 C.617€5007D—1é 
0.32188554D0-Cl 0.46862175D-Cé 
0.3545&815D-0€ 

0.20503732D Cl —-C.~17€1€2€8D-15 
0.118698C10-C1 0.212559140-Cé 
0.12640418D-0€ 

0.22591483D Cl 0.2644€274D-LE 
0.40572610D-C2 C.890S3523D-07 
0.1521814ED~G7 
0.24456866D O01 C.7241C0S41D-1€ 


Oseclicele2ol-0c 
= O.1>ECCCOCCU=05 


ACCURACY CESIRED HAS AQT 6EEN ACHIEVED, 
ITERATICAN CCATIANUES.. 


fede 





PR(xX=C),ECLNC 
Bote ERRCR 
XC,FXO,IFLAC 
PR(X=1),ECUNC 
EST. ERKRCR 
XC,sFXC,IFLAC 
PR(X=2),ECUNL 
EST. ERRCR 
XC,FXC,ITFLAC 
PR(X=2),@CUNC 
feetie ERRCR 
XCyFXC,IFLAC 
FR(X=4),ECUNC 
EST. ERKRCR 
XC»FXC,TFLAC 
PR(X=5),ECLNC 
EST. EFRCR 
XCsFXC,IFLAC 
PR(xX=€),2CUNC 
EST. ERRCR 
XCsFXO,IFLAC 


ACCURACY CESIREC ACHIEVEC AT 


TABLE XVIII 


ITERATICN # Ll 


GesoSsz2990eR CC 


C.95894C18D CO 
One s54652D C0 


O.12781526D Cl 
Orel eooro? DCC 


0O-15€34/776D Cl 
0./95328928D-Cl 


0.18195852D Cl 
GO .232188563D-C1l 


0.20503833D Cl 
0.118658C03D-Cl 


0.22591594D Cl 
0 .405/26040-Q2 


0.24486524D Q1 


he 


11TH 


Ge2742Z21(260-C6 


SU. bous2ceoap-O06 
~Os74 265 1540-1 7 


0.29224EEI(D-OE 


SOelaterso()—O¢ 


Gr CGC eo 12 
O.21C25E26D-0€ 
O.LC1S2154D-0€ 
Ga 055 cl 50-16 
Orcs (24CSU= CE 
0.94SE1C00D-C? 
Capt 241Ci 0-16 
Gsezsié2iCp-07 
0.20€2/1421D0-07 


Sie Hee cecal — 2 


0.28425182D-O07 


-0.81489224D-08 


C.1477S243D-1¢E 
CelisSel72S2)-Cil 


TOsioecelog0—-CE 
weno ecU 4a 15 


ITERATION 
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